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ABSTRACT 

We show that surface waves along interstellar current sheets closely aligned with 
the line of sight lead to pulsar scintillation properties consistent with those observed. 
By contrast with previously considered scintillation drivers, our mechanism naturally 
produces the length and density scales of the ISM scattering lenses that are required 
to explain the magnitude and dynamical spectrum of the scintillations. In our picture, 
the parts of warm ionized interstellar medium that are responsible for the scintilla- 
tions are relatively quiescent, with scintillation and scattering resulting from weak 
waves propagating along magnetic domain boundary current sheets, which are both 
expected from helicity conservation and have been observed in numerical simulations. 
The model quantitatively predicts the spacing and amplitudes of inverted parabolic 
arcs seen in Fourier-transformed dynamical spectra of strongly scintillating pulsars. 
Multi-frequency, multi-epoch low frequency VLBI observations can quantitatively test 
this picture. If successful, in addition to mapping the ISM, this opens the door to pre- 
cise nanoarcsecond pulsar astrometry, distance measurements, and emission studies 
using these lOAU interferometers in the sky. 
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1 INTRODUCTION: OBSERVATIONS AND 
THEORETICAL CHALLENGES 

Compact radio sources provide a precision probe of the ion- 
ized interstellar medium (IISM). The propagation speed of 
radio waves depends on the density of free electrons, and 
therefore the spatial inhomogeneity of the IISM may result 
in refractive and diffractive abberation and scattering. This 
causes s cintillation (time-variability) of the compact radio- 
sources (|Scheueilll963 : iBlandford et al.lll986l : iRickett et al.l 
l2006l ). Observations of the pulsar scintillations are partic- 
ularly interesting for infering the IISM properties, due to 
the brightness of many pulsars that have been observed for 
other purposes over long time intervals. 

However, despite of considerable effort, the small-scale 
structure of the IISM has remained enigmatic. At a first 
glance, one expects the scattering to be caused by den- 
sity inhomogeneities produced from turbulent motions of the 
IISM. However, this picture is in contradiction with the last 
decade of the pulsar scintillation data. There, a major obser- 
vational progress of th e ISM s c intilla tions has been achieved 
through the Stincbrin g et al.l (|200ll ) detection of parabolic 
structures in the Fourier-transformed dynamical spectra of 
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strongly scintillating pulsars. These parabolic structures im- 
ply that for these pulsars the radio-wave scat tering occurs 
most l y within one or se v eral thin screens (Walker et al.l 
l2004l : ICordes et all l2006l : IWalker et all 12008 ). Moreover, 
the multiple "inverted parabolae" I Hill et al. |200d) show 
that the scattering inside the screen is strongly inhomo- 
geneous and occurs i n localized clumps. The latter was 
recently confirmed by iBrisken et all l|2010l) who have ob- 
tained the VLBI scattering image PSR B0834-I-06. They 
have found that not only the scattering image was clumpy 
but that the clumps lined up along a thin line. All of 
these observational facts present a major challenge for the- 
oretical interpretation where the scattering is caused by 
densi ty inhomo geneities produced in a turbulent cascade 
(Lith wick fc Go ldrcich 2 0011 ). In particular, (1) the origin 
of the screens is unexplained, and (2) the strongly non- 
gaussian scattering requires AU-size regions which are over- 
press urized by factors of ^ 10^ relative to the ambient warm 
ISM (Romani et al. 1987"); no conventional physical mecha- 
nism has been proposed for how such regions may be formed. 

In this paper, we advocate a scenario in which the scat- 
tering is produced by a refractive structure, with turbulence 
playing a secondary role. Radio- wave scattering by non- 
turbulent large-sca le refractive structur es has been previ- 
ously considered bv lRomani et al.l l|l987i ). mostly in the con- 
text of the so-called extreme-scattering events observed by 
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iFiedler et all (|l987h . Recently, iGoldreich fc Sridhaj (|2006l ) 
proposed that the image of the SgrA* radio-source is 
strongly scattered by several reconnection sheets that are 
closely aligned with the line of sight to SgrA*. In this pa- 
per, we deve lop further the Go ldreich & Sridhar's (2006) 
idea [see also lPen fc Kinel (|2012l )] and apply it to construct 
a quantitative picture of pulsar scintillations. Namely, we 
consider a scenario where the pulsar radio-wave scattering 
occurs due to several weakly corrugated reconnection sheets 
that are closely aligned with the line of sight to the pulsar; 
see figure [T] We show that this scenario provides compelling 
explanations for previously unexplained features of the scin- 
tillations: (1) the "scattering screens" are simply effective 
descriptions of such sheets; their location is marked approx- 
imately by the sheets' intersections with the line of sight, (2) 
"the scattering clumps" correspond to those parts of sheet 
folds where the sheet is parallel to the line of sight; the 
strength of the scattering follows a strongly non-Gaussian 
distribution, even though the corrugation itself is assumed 
to be a realization of a Gaussian distribution, and (3) the 
strong non-isotropy of the Brisken et al. 2010 image is a con- 
sequence of the sheet's inclination, with the clump locations 
aligned along the direction perpendicular to sheet's line of 
nodes. The plan of the paper is as follows: in the next section 
we briefiy describe the origin of the reconnection sheets, in 
section 3 we derive the model for the fold statistics, in sec- 
tion 4 we derive the lensing by the corrugated sheet, and 
in section 5 we compute the Fourier-transformed dynami- 
cal spectrum and demonstrate the parabolic structures. In 
section 6 we conclude. 



2 ASTROPHYSICAL PICTURE 

2.1 Two Regimes of Lensing: Diffractive vs 
Refractive 

Two regimes to generate pulsar scintillation have been con- 
sidered. In the diffractive regime, the scattering/lensing an- 
gle is determined by the ratio of the wavelength A of the 
radio-waves, and the characteristic spatial scale D of the 
scattering structure: ~ X/D. The brightness of the scat- 
tered image is determined by the amplitude of the wavefront 
modulation caused by the scattering structure. To explain 
the observed angles in the range 1 — 100 mas at wavelength 
~ Im, requires structures in the ISM on scales of order 
l(f~^m. This imposes unexpected properties on the IISM, 
since this scale much smaller than the coloumb mean free 
path of protons and thus compressive perturbations of this 
scale are overdamped and decay esponentially on the sound- 
crossing time-scale. If one still assumes that somehow these 
perturbations are created and maintained, then the angular 
image of a pulsar, and therefore its dynamic spectrum, are 
superpositions of thousands of AU/10*m weak structures 
(for a 1-d scattering image), and are expected to be roughly 
Gaussian. The length scale is given by the path length differ- 
ence, which in turn is inferred from the inverse scintillation 
correlation frequency. The asymmetric parabolically struc- 
tured 2-D power spectrum of the dynamic spectrum, and 
the VLBI image of the scattering disk consisting of several 
prominent clumps, are inconsistent with such a picture. The 
number of 10^~*m eddies along the line of sight can be large, 



possibly 10^^ for a pulsar at a distance of ~ kpc. The total 
scattered power comes from the cumulative projected vari- 
ation in refractive index, which grows as the square root 
of the number of eddies. Each eddie only needs to change 
the refractive index by a part in a million of the cumulative 
change, so a refractive index variation of a part in ~ 10^^ 
could account for the strong scintillation. But the super- 
position of scattering from such a large number of eddies 
would surely look very Gaussian by the central limit theo- 
rem. A requirement for the sum of 10^'^ contributions to ap- 
pear intermittent in projection would lead to unphysically 
overpressurized eddies. 

A second mechanism is due to refractive lensing. In this 
scenario, the bending angle is determined by Snell's law, i.e. 
the change in refractive index and the angles of incidence. 
The refractive index of a plasma is n = 1/ -^/l — ujp/to^, with 

plasma frequency ujp = y^Uee^ /eoirie- Pulsar observations 
are done at frequencies much higher than the plasma fre- 

quency, for which we expand n — 1 ~ = 1.8 x 10~*ne at 
wavelengths of a meter. The observed scattered images at 20 
mas require deflection angles of at least 40 mas, correspond- 
ing to He ~ 10/cm^ neglecting geometric alignment factors. 
However, the mean density of the IISM is determined to be 
of the or der 10~'^, as measured fr om the pulsars' dispersion 
measure dLorimer fc Kramei|[200i ). Therefore, the refractive 
picture is also challenging to reconcile with the data, since 
the observed scattering angles would naiVely require frac- 
tional changes in free electron density of ~ 10'' which are 
difficult to understand or confine. This contraint, however, 
is alleviated when one considers refra ctive sheets that are 
close ly aligned with the line of sight (|Goldreich fc Sridharl 
l2006h . 

Historically, refractive lensing was used to interpret long 
time variability, and diffraction for the minute time scale 
effects. Recently, it was understood that the refractive im- 
ages result in an interference fringe pattern (I Walker et al.l 
2004) which have similar time and frequency scales for flux 
modulation as diffractive effects. The frequency and time 
scaling was historically interpreted as related to an under- 
lying stochastic diff r active process driven by turbulence. 
IGoldreich fc Sridhail ()2006l ) showed that refractive lensing 
by aligned sheets results in scintillation similar to the diffrac- 
tive picture. 



2.2 Physical origin of the reconnection sheets 

The interstellar medium is stirred on scales of parsecs by 
various energetic processes, including supernovae, ioniza- 
tion fronts, spiral density waves, and other phenomena. 
These processes are generally short lived, and we conjec- 
ture after the stirring, the warm medium relaxes to an near- 
equilibrium configuration on small (several AU) scales. Cur- 
rent numerical simulations of the supernova-driven turbu- 
lence in the warm ISM of the Galaxy (e.g.. Hill et al. 2011) 
do not have the resolution to tell how realistic this assump- 
tion is. In the presence of helicity, the equilibrium magnetic 
fields are confi gured a s interlaced twisted tori, which are 
lo ng lived (IBraithwaite fc Spruit 2004). It has been shown 
bv lGruzinovl ( 20091 ) that a "generic magnetic equilibrium of 
an ideally conducting fluid contains a volume-fllling set of 
singular current layers." In this picture, the magnetic fields 
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are locally almost parallel, with discontinuous interface re- 
gions, a bit like magnetic domains in a ferromagnet. Singu- 
lar current sheets have al so been seen in the simulations of 
ISchekochihin et al.l (|2004l ). 

At the boundary between between magnetic field con- 
figurations, current sheets maintain the discontinuities. De- 
pending on the nature of reconnection, these current sheet 
may be self-enforcing due to inflow of fresh fields, maintain- 
ing a thickness potentially as thin as an electron gyromag- 
netic radius. At constant temperature, the pressure equilib- 
rium in the direction perpendicular to the sheet implies that 
the electron density inside the sheet is enhanced by a factor 
R given by 



R- 



Bl 



4^Pout ' 

where Bin and Brmout is the magnetic field in and outside 
the sheet, and Pout is the non-magnetic pressure outside the 
sheet. The right-hand side of the above equation is thought 
to be of order 1 for the IISM, but may be significantly larger 
if the IISM is magnetic-pressure dominated. 

Depending on the ratio of heating time due to reconnec- 
tion to cooling time, the thin sheet may have an enhanced 
temperature, in which case it can be underdense, with 7? — >■ 
in the limit of a strong entropy injection. 

The lensing geometry is shown in Figure [1] A crucial 
ingredient in our model is that the reconnection sheet is 
assumed to be weakly corrugated. Let is assume that the 
corrugation pattern is fixed, and vary the inclination angle. 
In this case, as can be seen from figure[T] there is character- 
istic inclination angle a at which the perturbations in the 
sheet can generate caustics in the projected surface density. 
This angle a is the ratio of the wave peak displacement to 
its wavelength. We envision ranges of a ~ 10"=^ - 10"^ 
When the caustics are present, they become effective refrac- 
tive scattering centers for the radio waves. The resulting 
lense features multiple images that closely line up along the 
direction perpendicular to the line of nodes of the scattering 
sheet. 



2.3 Surface Dynamics 

We assume that the current sheet is physically thin, < AU. 
Theoretically, the thickness of current sheets is not under- 
stood, so we choose this scale to be thin enough to explain 
the smallest scale observed structures. On each side the mag- 
netic field points in a different direction. T he change in al- 
venic properties gives rise to surface waves (|jain fc Roberts! 
Il991 : Joardcr ct al. 2009), whose amplitude decays exponen- 
tially with the distance to the current sheet and which are 
mathematically analogous to deep water ocean waves. The 
restoring force is due to the difference in magnetic field com- 
ponent projected along the wave vector. Like ocean waves, 
these waves penetrate about a wavelength into each side. We 
will be considering wavelengths of hundreds to thousands of 
AU whose projected wavelength is related to the observed 
lensing structures ~several AU, so the thickness of the cur- 
rent sheet itself is neglible as far as the dynamics of the 
waves are concerned. Seen in projection along the aligned 
sheet, the projected wavelengths will be ~ AU, reduced by 
the alignment angle a. 

The displacements are transverse to the wave vector. 



and perpendicular to the sheet. While alvenic in nature, 
the surface modes possess only one polarization, unlike bulk 
waves. We speculate that this allows such modes to be long 
lived, since the single polarization natu re protects them from 
the n ormal MHD turbulent cascade (jGoldreich fc Sridharl 
Il99'it l. These waves resemble a flag blowing in the wind. 
Disturbances travelling along the sheet are decoupled from 
bulk waves. Being confined to a sheet, the amplitude away 
from a source drops as oc 1/r instead of the normal inverse 
square law. The 2-D analogy to Olber's Paradox leads to 
flux being dominated by the cumulative efi'ect of far away 
sources, like waves on an ocean beach. A amplitude of order 
oAwavG, or about 10~'^-10~^ of the wavelength of the sur- 
face wave, is sufficient to cause the sheet to appear folded in 
projection; here Awavo is the wavelength of the surface wave. 



3 FOLD STATISTICS 

We expect a minimum wavelength of these surface waves, 
which is substantially larger than the thickness of the sheet. 
Short wavelength perturbations are not bound to the sur- 
face, and can dissipate into the bulk. The exact cutoff de- 
pends on unknown factors, including the distance to the 
source, and vertical structure of the current sheet. Non- 
linear effects also cause short wavelength waves to dissipate, 
just like sound waves in the air. Primarily waves with am- 
plitude larger than a contribute to scattering. The strength 
of damping depends on the distance to the source in units of 
wavelength, with shorter wavelengths being more damped. 

We model the waves as a displacement function ({x) 
which is a Gaussian random field with a correlation function 
that is a Gaussian, 



(2) 



C(r) = (C(x)C(a; + r)) = exp - 

where A is the mean amplitude of the displacement and a is 
the surface-wave coherence scale. We denote the projected 
coherence scale (Jy = aa. Figure [5] shows a realization of a 
sheet with a random fiuctuations. We used an inclination 
slope a = 1/8, and a correlation length along the sheet of 
350 units, which is ay ~ 44 grid units in projection. The 
displacement amplitude is chosen as j4 = 40. The dimen- 
sionless fluctuation 5 = A/a ^ 1 is about unity, meaning a 
1 — a flucuation can result in a caustic fold, which we will 
discuss further below. 

The current sheet corresponds to a change in magnetic 
and thermal pressure and thus has a refractive index differ- 
ent from the ambient ISM. For simplicity, we assume here 
that the sheet has constant thickness, and consider its opti- 
cal depth as is relevant for refractive lensing. In projection 
along the line of sight, the column density of the sheet re- 
sults in a highly non- Gaussian distribution. Figure |3] shows 
the column density distribution in a simulation. Folds oc- 
cur when the gradient of the displacement is equal to —a, 
which as described above is quantified by the fluctuation 
amplitude 5. The value chosen here makes folds common, 
occuring roughly once per correlation length, and yet make 
multiple folds overlapping in projection rare. Each fold re- 
sults in two caustics, with characteric separation ay. This 
is wel l understood from the theory of extrema of surface 
waves l|Longuet-Higginsl Il957l ) . 
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Figure 1. lensing geometry. The earth is at left, pulsar at right. A section of the scattering sheet is in the middle. The dashed line shows 
the unperturbed light path. The dotted line shows the path of an image lensed by a fold caustic of the projected sheet. 
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Figure 2. Sheet with transverse perturbations. The upper panel 
shows an zoomed version of a short section. 



As before, a is the angle between the screen and the 
line of sight. Then the optical depth is p oc 1/a for a <^ 1, 
and Pscreen(p) OC 1 /p^ , where Pscreen{p) is the probability of 
a piece of screen to have the optical depth of p. However, we 
are interested in the probability density with respect to the 
impact parameter relative to the line of sight, and not with 
respect to the location on the screen. For nearly aligned 
screens, this is not the same thing. In particular, part of 
the screen with low a occupies less of the impact-parameter 
space than the part of the screen with the same area but 
high a. Thus, 



-Pimpactparamctcr (p) 1/p 



(3) 



which is a strongly non-Gaussian distribution. 

Equation ([3]) describes the 1 point PDF in figure [S] The 
deflection angle is determined by the gradient of the density, 
P'. 

As in iPen fc Kind l|2012D . we use the notation of grav- 
itational lensing. The phase delay through the lens is de- 
scribed by the lensing potential 



Uedze 



(4) 



where the effective density = p/dz. The lensing depends 
only on the product of density and thickness. The conver- 
gence K = dgTp/2 is given by the second angular derivative 
of the potential, and thereby the projected density p. The 
magnification gives the brightness of the image. In 1-D it is 
yit = 1/(1 — 2k), derived from the determinant of the ampli- 
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Figure 3. projected density. The upper panel is a zoom of the 
central portion of the lower panel. Each time the sheet folds in 
projection, we see a double peaked caustic structure in density. 
The characteristic separation between the two peaks is a projec- 
tion correlation length cry, in this case 44 grid units. 



fication matrix. A negative amplification corresponds to a 
flipped image of odd parity. 

The regions near the locations where the screen is par- 
allel to the line-of-sight, which the call the caustics, give 
rise to the localized clumps in the pulsar's scattering im- 
age, which produce the inverted parabolic arcs in pulsar 
secondary spectra. 



4 LENSING 

T he lensing of this den sity sheet can be computed in analogy 
to lPen fc Kin3 (|2012l ). 

Given the projected column density distribution in Fig- 
ure |31 and assuming that the sheet has a fixed optical depth 
along its normal, we can compute the mapping of apparent 
angle on the sky to position in the source (pulsar) plane. 
The caustics in the projected density distribution lead to 
large angle deflections, and multiple images, whenever the 
sheets are aligned closely enough for the caustics to form. 
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Figure 4. deflection angle mapping. The horizontal axis is angle 
on the sky, and the vertical axis is the intersection of this light 
ray on the source plane. Whenever multiple different directions 
on the sky intersect on the same position in the source plane, 
multiple images are formed, which form a coherent interference 
pattern. 



This explains why only a small fraction of the current sheets 
contribute to scintillation. 

The system depends on the dimensionless parameter 5, 
the ratio r of sheet thickness z to the projected correlation 
length ay, and the ratio of index of refraction change to the 
angular size of the projection correlation length. 

The model predicts the number density of images and 
their fluxes as a function of their angular separation from 
the Une of sight (and therefore as the function of time lag). 
It has a small number of dimensionless parameters: S, r, C 
(defined below). 

There is a dimensional scaling of time units, which is a 
function of pulsar transverse velocity, projected screen size, 
and observing frequency. This is generally parameterized as 
the DISS time scale, tDiss ~ {^/D){L/v) where D is the 
size of the lensing region, L is the distance to the pulsar, A 
is the observing wavelength, and v is the pulsar transverse 
velocity. Note that in this model, the lensing is refractive, 
sharing the time scales from diffractive models, but none of 
the physics. 

We show a histogram of image magnifications in Fig- 
ure [51 which can be compa red to holographic flux measure- 
ments I Walker et al.l [20081 ). The lensed images correspond 
to the inverted arclets seen in secondary spectra. The model 
predicts the number of images as a function of separation, 
shown in figure [S] 

The separation of images is related to the separation of 
caustics. For the "common" limit of our simulation, 5 ~ 1, 
this is given by the projection correlation length ay. For 
5 <C 1, the angular density of caustics becomes very rare, 
proportionate to an error function. 

The number of images is determined by the number of 
light folds along the line of sight, i.e. how often the deflection 
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Figure 5. PDF of image magnifications. The peak near 1 are 
images at the unscattered positions, while the population on the 
left are lensed images. The peak occurs at roughly 
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Figure 6. PDF of image positions per logarithmic distance in- 
terval on the sky y. The projected correlation length is <jy ~ 44. 
The characteristic convergence k ~ 20, leading to a cutoff near 
UyK ~ 1000. 



angle is larger than the separation to the line of sight in 
Figure [4] 

The max;imum deflection angle of an image is the change 
of refractive index in the sheet, amplified by the alignment 
angle on a fold caustic. For a sheet much thinner than 
the projected fluctation amplitude, the maximum projected 
density enhancement C = p/p relative to the projected mean 
sheet density p ~ riedz/a at a caustic is the square root 
of the radius of curvature to the thickness of the sheet t, 
C = ayl\fAT in the limit a"^ j A » r. 
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Figure 9. secondary pulsar spectrum. The inverted parabolic 
arcs arise naturally in this model. 



images forms the inverted para bolic arcs, qualita t ively simi- 
lar to those that are observed in lStinebring et al.l (|200ll ) and 
iHill et all (|2005l ). 



Figure 7. Average Spatial Distribution of flux per logarithmic 
distance interval. In a Gaussian model, the flux per image drops 
as 1/0. 




Figure 8. dynamic pulsar spectrum. Horizontal axis is time, ver- 
tical axis is frequency. We reproduce the characteristic criss-cross 
pattern observed in real scintillation spectra. 

5 SIMULATED DYNAMIC SPECTRA 

With the density field, we can solve the lens equations to 
simulate dynamic spectra. By adding the voltages on each 
image with their appropriate amplitude and phases, we sim- 
ulate the dynamic spectrum, shown in figure [S] 

A 2-D fourier transform maps this dynamic spectrum 
into a secondary spectrum, shown in figure |9] 

We find that the interference of these discrete, co- linear 



6 DISCUSSION 

We can estimate the length scales involved in making the 
current sheet that would produce the observed scintillation 
pattern. This theory requires as input a current sheet thick- 
ness, inclination angle, curvature, amplitude of waves, and 
dissipation scale. 

The thickness of the sheets can be estimat ed by consid- 
ering the magnification of images. As shown in I Pen fc Kinj 
(2012), the fiux is roughly the thickness divided by the im- 
pact parameter. This follows from flux conservation of lens- 
ing: the net flux is conserved, and flux changes by order 
unity at impact parameters of order the physical size of the 
lens, so the typical fiux off-axis is roughly the ratio of the 
furthest distance at which at image forms, to the size the 
lens. The caustic itself only contains a small fraction of the 
integrated flux. 

The projected wavelengths are observed to be of order ~ 
10 AU; this suggests a typical thickness of h < 0.1 AU to ex- 
p lain t he ~% scattering intensities observed in lBrisken et al.l 
( 20ldi) . 

The largest observed deflection angles at meter wave- 
length are 7 ~ 0.01", which requires an electron density 
change Sue ~ lOOa/C cm~^. For a typical interstellar 
plasma density of rie ~ 0.03 determined from pulsar disper- 
sion, we need an alignment of a fraction of a degree, with 
fiuctuation amplitude C ~ 30 and a wavelength of A ~ 1000 
AU. This combinat ion is not un i que. 

As discussed in lPen fc Kind (|2012l ). the phenomenology 
of the Extreme Scattering Events prefers underdense lenses. 
In an underdense sheet, the maximal change of density is 
the density itself. With the assumption Srie ~ rie , we obtain 
a ~ 10~^. The probability of seeing a sheet at such an angle 
is ~ a^, requiring the existence of ~ sheets along the 

line of sight: if we imagine each "sheet" to be a thin disk, the 
probability of seeing a disk at angle a gets one contribution 
from the intrinsic alignment distribution, and one more from 
the reduced geometric cross section of an aligned sheet. 

For pulsar B0834-I-06, the distance is ~ 0.64 kpc (as de- 
termined from the Dispersi on Measure, and con sistent with 
VLBI and ISM geometries l|Brisken et alll2010l '). also direct 
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distance of order of the effective emission height multiphed 
by the ratio of the pulse width to the rotation period. Using 
* VLBI mapping of the scattering geometry, one can precisely 

predict the change of scintillation pattern as a function of 
\\\ III 11 \' I pulse phase. The effective astrometric precision can be sub 

nano arcsecond. 



Figure 10. Simulated scattering image of the pulsar. The pul- 
sar is modeled as a unit disk, with a much exagerated scale. Flux 
conservation requires that the sum of image areas equals the orig- 
inal disk. Th i s com pares favourably with the VLBI images of 
iBrisken etahl 1 I2OIOI) . 

parallax: Deller and Brisken, unpublished ), giving a typical 
sheet separation of s ~ 0.1 pc. 

These estimates are qualitative. One expects current 
sheets to come in a range of sizes, curvature and perturba- 
tion amplitude. The thickness might also vary. 

One of the attractive features of our mechanism is 
that it explains very naturally the 1-d image of Brisken 
et al. (2010). The reduced dimensionality of the scattering 
image comes from the fact that the deflection created by 
the screen is mostly in the direction perpendicular to the 
screen's line of nodes. Therefore, the scattering clumps will 
also form a line that is perpendicular to the line of nodes. 
This is demonstrated in Fig. 1101 where we show a simulated 
2-d scattering image of a pulsar. The alignment of the scat- 
tering clumps is apparent. 



7 FUTURE POTENTIAL 

In our picture, pulsar scintillation is dominated by a small 
number of magnetic discontinuities highly aligned to the line 
of sight. Surface waves will move very slowly in this pro- 
jected geometry, allowing for a precise determination of the 
geometric properties. This allows the use of these sheets as 
lenses to study both pulsars and the ISM (|Pen et al.|[2013l ). 

7.1 ISM dynamics 

Free parameters in our model include the thickness of the 
sheet and the inclination angle. These can be inferred by 
broad band measurements of pulsar scintillation, as follows. 
Firstly, the apparent position of images is expected to shift 
by a distance of order of the sheet width, as one decreases 
the observing frequency from the highest critical frequency 
at which the image forms, to a factor of two below. Secondly, 
the co-linearity of the images is related to the inclination 
angle of the sheet: the more aligned is the sheet with the 
line of sight, the greater is the aspect ratio of the scattering 
image. 

7.2 Pulsar Emission imaging 

A straightforward application is the study of the reflex mo- 
tion of the emission region of the pulsar across the pulse 
phase. One expects the apparent emission region to move by 



7.3 Distance Measurement 

It is tempting to use these lenses to obtain precision paral- 
lax distances to pulsars. The challenges is to keep the lens 
stable over a year, when typical pulsars move by much more 
than an astronomical unit, making the differential measure- 
ments challenging. In the sheet lensing scenario, one can 
imagine obtaining widely separated scattering images: the 
lensing angle scales as oc A^, so at low frequencies, for ex- 
ample with LOFAR LBA, hundreds of AU are probed. Over 
the course of a year, persistent images should remain visible, 
and the interference pattern could be studied as a function of 
annual modulation. This would result in direct parallax dis- 
tances with nano arc second precision, enough to determine 
pulsar distances for coherent gravitational wave detection 
(,Bovle fc Penii2012>'l . 

8 CONCLUSIONS 

We have presented a quantitative theory of pulsar scintil- 
lation inverse parabolic arc s. T hese extend r e cent i deas of 
iGoldreich fc Sridhad (|2006l ) and |Pen fc Kind l|2012l ) about 
thin current sheets as the scattering objects in the ISM, 
which naturally explain the large angle scattering observed 
in pulsars and some extragalactic sources. 

This picture could explain all scintillation phenomena 
from refractive lensing, all with structures greater than 0.1 
AU. The apparent diffractive structure results from the in- 
terference between refractive images, and no diffractive scat- 
tering is needed. 

In this scenario, VLBI monitoring of pulsars on time 
scales of weeks, at multiple low frequencies, can allow fore- 
casts of the scattering behaviour. This in turn could improve 
gravitational wave timing residuals. The same scenario also 
enables the coherent use of the scattered images as a gigantic 
interstellar interferometer to map the motions of the pulsar 
emission regions. 
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